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We employ the exact eigenstate basis formalism to study electrical conductivity in graphene, in 
the presence of short-range diagonal disorder and inter-valley scattering. We find that for disorder 
strength, W > 5, the density of states is flat. We, then, make connection, using the MRG approach, 
with the work of Abrahams et al. and find a very good agreement for disorder strength, W = 5. 
For low disorder strength, W = 2, we plot the energy-resolved current matrix elements squared 
for different locations of the Fermi energy from the band centre. We find that the states close to 
the band centre are more extended and falls of nearly as 1/Ef as we move away from the band 
centre. Further studies of current matrix elements versus disorder strength suggests a cross-over 
from weakly localized to a very weakly localized system. We calculate conductivity using Kubo 
Greenwood formula and show that, for low disorder strength, conductivity is in a good qualitative 
agreement with the experiments, even for the on-site disorder. The intensity plots of the eigenstates 
also reveal clear signatures of puddle formation for very small carrier concentration. We also make 
comparison with square lattice and find that graphene is more easily localized when subject to 
disorder. 



I. INTRODUCTION 

Graphene has been long studied as a theoretical toy 
model not only to understand it's appealing physical 
propertieSfi"— but also as a basic building block of vari- 
ous Carbon allotropes like graphite^i^ and more recently 
fullerenes and nanotubesii While graphite is the three 
dimensional allotrope of Carbon and could be formed 
by the Bernal stacking of graphene sheets, fullerene and 
nanotubes are the zero and one dimensional allotropes, 
formed by introducing pentagonal impurities and rolling 
the graphene sheets, respectively. After its experimental 
isolation in 2004;^"— there has been a renewed interest in 
studying various properties of graphene sheet, both the- 
oretically and experimentally, as well as due to potential 
technological applicationsiiiii^ 

Graphene consists of a single sheet of Carbon atoms 
arranged on a honeycomb lattice. Basic properties of 
graphene are well described by a tight-binding model for 
the TT-orbitals which are perpendicular to the graphene 
plane at each Carbon atom. The effective low-energy the- 
ory states that the charge carriers in graphene are mass- 
less Dirac fermions, characterized by a linear dispersion 
relation and a linear energy dependence of the density 
of states which vanishes at the Fermi level implying a 
semi-metallic behaviour for graphenci ^'^'^'^ Graphene has 
attracted a lot of attention recently not only due to its 
potential technological applications but also for under- 
standing of physics in 2D systemsiii^. Its low energy de- 
scription mimics (2+l)-dimensional quantum electrody- 
namics and hence graphene could act as a testing ground 
for various relativistic phenomenaJ^ 

Early experiments on graphene have revealed that the 
conductivity at low temperatures is directly proportional 
to the carrier concentration (or gate voltage) except for 
very low carrier concentration. For zero gate voltage, 
the conductivity approaches a robust minimum univer- 
sal value proportional to e^/hM^ This could not be ex- 
plained by the Born approximation which predicts a con- 



ductivity independent of carrier concentrationiii^ Other 
interesting properties include anomalous integral quan- 
tum Hall effect and suppression of weak localization.— 
Recent experiments, however, show that the depen- 
dence of conductivity on carrier concentration could 
vary from sub linear to superlinear for different carrier 
concentrations . It has been argued that presence of 
impurities in graphene is the main contributor towards 
its electronic properties'^ The importance of disorder in 
graphene could most easily be emphasized by observing 
that the universal conductivity suggested by the theo- 
retical studies on defectless graphene sheet is 2-20 times 
smaller than the observed conductivity close to the Dirac 
pointsjii 

The Boltzmann conductivity for graphene is given 
by ctq = {e^ /h){2EpTi^/h). The observed conductivity 
rises linearly with carrier concentration in graphene and 
N{Ep) (X y/n , where N{Ep) is the density of states at 
the Fermi energy and n is the carrier density. This im- 
plies that the scattering rate, T{Ep) oc ^/n. On the other 
hand, for weak local scatterers. Born approximation pre- 
dicts r"-'^ cx ninipN{Ep) oc ^/n where nimp is the im- 
purity concentrationjii In view of this discrepancy, vari- 
ous investigations, both theoretical and numerical, have 
been carried out in order to understand the behavior of 
graphene under various types of disorder j^i-'^^'^^^'^^i'^^i^"^" — 
such as vacancies,—!^ charged carriers j^^^"— on-site 
disorder;^!^ long range on-site disorder;^ off-diagonal 
disorder off-diagonal disorder with sign change prob- 
ability in the hopping termi^ 

Vacancies have been proposed to induce localized 
states, extended over many lattice sites, which are sen- 
sitive to the electron-hole symmetry breaking^^ Detailed 
studies in the presence of both compensated and uncom- 
pensated defects reveal that they could modify the low 
energy spectrum in graphene drastically like there could 
be quasi-localized zero modes and introduction of gap in 
the DOS^ 

For charged scatterers, Nomura et ali^ have argued 
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on the basis of Boltzmann transport theory that the hn- 
ear dependence of conductivity on carrier concentration 
could be explained. They find that states close to the 
Dirac point are delocalized leading to cr^i^. Also, one 
could observe antilocalization if the inter-valley scatter- 
ing is weak. On the other hand, if inter-vally scatter- 
ing is large, all states could be localized due to accu- 
mulation of Berry phases. Conductivity in the presence 
of random charged impurity is also studied by Hwang 
et alM' They find linear dependence of conductivity on 
carrier concentration for high carrier density. However, 
for low carrier density, they argue, that system devel- 
ops some inhomogenities (random electron-hole puddles) 
which implies that this domain is dominated by localiza- 
tion physics. They also conclude that change of bias volt- 
age may change the average distance between graphene 
sheet and the impurity in the substrate which could 
lead to sub- and super-linear conductivity dependence 
on carrier concentraion. Pereira et al^ have argued that 
there could be a "critical coupling" distinguishing strong 
and weak coupling regimes in the presence of unscreened 
Coulomb charges. They also find bound states and strong 
renormalization of Van Hove singularities in the DOS. 
Vyurkov et a/i^have argued that the intrinsic conduc- 
tivity of graphene (ambipolar system) is dominated by 
strong electron-hole scattering. It has a universal value 
independent of temperature. It is shown that conduc- 
tivity could be proportional to V or V'^ depending on 
the other scattering mechanisms present like those on 
phonons by charged defects. In the unipolar system, it 
is argued that electron-hole scattering is not important 
and conductivity is proportional to V^^""^. Trushin et al^ 
have reported an analytical calculation for Boltzmann 
conductivity with screened Coulomb scatterers and both 
electron- hole coherent and incoherent solutions. They 
find that the experimentally observed dependence of con- 
ductivity on n could be explained by the electron-hole 
coherent solution. 

For diagonal disorder, it is found that beyond a cer- 
tain threshold disorder strength, bound states appear 
beyond the band continuum and resonant states could 
appear at low energies. In the infinite disorder strength 
limit, results match with that of vacancies. In the off- 
diagonal disorder case, strong low-energy resonances ap- 
pear. There are, however, no bound statesi^l Localiza- 
tion studies on graphene with on-site disorder carried out 
by Xiong et al^ reveal that all states are localized in the 
case of random diagonal disorder which is consistent with 
Anderson localization. Lherbier et al^ have also stud- 
ied the energy dependent elastic mean free path, charge 
mobilities and semi-classical conductivity in the pres- 
ence of Anderson-type disorder, using real space order 
N Kubo formalism, for both two-dimensional graphene 
and graphene nano-ribbons (GNRs). It was found that 
the systems undergo a conventioanl Anderson transition 
in the zero temperature limit. 

Lewenkopf et alM- have studied long range diago- 
nal disorder (Gaussian scatteres) at finite concentration. 



They find that conductivity increases as disorder strength 
is decreased. It is shown that conductivity depends only 
on disorder strength and ratio of the system size to dis- 
order correlation length. This dependence could vary 
between sublinear to superlinear depending on disorder 
strength. For fixed disorder strength, conductivity in- 
creases with doping concentration. In the presence of 
strong long-range impurities, Zhang et a/j^ have shown 
that states close to the Dirac points are localized for suffi- 
ciently strong disorder strength and Kosterlitz-Thouless 
transition between localized and delocalized states is pro- 
posed which is seen in terms of the "current flow vector". 

In the case of random off-diagonal disorder (hopping 
disorder), localization studies by Xiong et al^ reveal 
that states close to the Dirac point are delocalized due 
to chiral symmetry. They find that the off-diagonal dis- 
order leads to a shape-dependent conductivity depending 
on the length to width ratio. However, if a sign change 
probability is introduced, they find that the conductivity 
becomes shape independent. Amini et al^ have shown 
that for disorder strength less than the hopping strength 
i.e. W < t, there is no localization. Disorder W slows 
down the Dirac quasi particles but preserves their na- 
ture. For W ^ t, localization sets in for states close 
to the Fermi energy, gap at energy close to the Dirac 
point and for W > t existence of mobility edge is pro- 
posed which starts at Fermi energy and moves towards 
the edges. States close to the Fermi energy are extended. 
They also propose the existence of disorder induced gap 
defined as the distance between the upper and lower mo- 
bility edges around the Fermi point. 

It should be noted that in all these works, the inter- 
valley scattering is assumed to be very small and hence 
not contributing towards conductivity. However, Klos 
et al^ have done a comparative study of conductivity 
using the tight-binding(TB) Landauer approach and on 
the basis of the Boltzmann theory and find a discrep- 
ancy between that results obtained by TB calculation 
and Boltzmann approach. 

Despite all the efforts, the issue of localization in 
graphene is currently highly debated from a theoretical 
standpoint. The obsered minimal conductivity (Tmin ~ 
Ae^ /h over a range of mobilities remains to be fully un- 
derstod. In view of a recent experiment by Ponomarenko 
et al^ and Katoch et al^ , which explores the dominant 
scatterers in graphene, and Horngei al^, which measures 
the high-frequency conductivity in graphene, unlike be- 
lieved so far, the primary reason for the linear rise of 
conductivity with carrier concentration is also debatable. 

In the present work, we will investigate the finite- 
frequency electrical conductivity and localization prop- 
erties of graphene in the presence of diagonal (on-site) 
disorder for various disorder strengths. Our exact- 
eigenstates approach implicitly takes into account the 
inter- valley as well as the intra- valley scatterings. This 
paper is organized into eight sections. 

In section II and HI, we introduce the two sub-lattice 
basis and evaluate the exact single-particle Green's func- 
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tion within the t-matrix approach for a single impurity on 
either sublattice. In section IV, we consider disordered 
graphene. Disorder is introduced via random fluctuation 
of the on-site energies of the 7r-orbitals. 

In section V, we use the Kubo-Greenwood formula to 
calculate the frequency-dependent conductivity for dif- 
ferent disorder strengths and for different system sizes. 
The system-size dependence is employed to perform a 
renormalization group analysis, and for moderate disor- 
der strength contact is made with the weak localization 
result of Abraham's et at. Later on, focus on the weak 
disorder strength. We study the frequency dependence 
of the averaged current matrix element squared for dif- 
ferent locations in the band (Fermi energies) and hence 
calculate conductivity and mobility for low-disordered 
graphene samples. Wc also study the dependence of en- 
ergy resolved current matrix elements squared on system 
size and disorder strength. 

In section VI, we study the average current matrix el- 
ements squared over a range of disorder strength and for 
different system sizes in order to gain a better insight 
into disorder induced localization in graphene. 

In section VII, we make a comparision between 
graphene and square lattice. We compare the normal- 
ized conductivity between the two. We also show the 
intensity plots for both the lattices for different values of 
disorder strength. 

Finally, our conclusions are presented in section VIII. 



II. FORMALISM 

As already mentioned above, graphene is a single layer 
of Carbon atoms on a honeycomb lattice. A honeycomb 
lattice, which is not a Bravais lattice, could be viewed as 
a composed of two kinds of sub-lattices. Nearest neighbor 
hopping takes electron from one sub-lattice to another. 
We start with the tight-binding(TB) Hamiltonian: 



(1) 



where the sum is over the nearest neightbours. This could 
be written in the two sub-lattice basis as: 



(5(k) 
S*{k) 

where, the hopping term in k-space is: 

S{k) = -t[e'''^^ + 2e(-^'=^'^/2)cos(^^^)] 
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It should be noted that this Hamiltonian is (2 x 2) and 
mixes the Dirac points. The energy eigenvalues are given 

by 
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FIG. 1: (color online) Density of States in graphene: (a) 
Density of States and (b) Real and Imaginary parts of 
the Green's function. 



In this basis, the low energy spectrum provides six Dirac 
points in the extended Brilluoin zone. These points are 
located at fc = (±27r/3, ±27r/3V3) and (0, ±47r/3V3). 
Expanding around the Dirac points, we find the linear 
dispersion for carriers in graphene. 

The free particle Green's function is given by the ex- 
pression: 



G°{k,uj)-- 



1 



'u}-H{k) 



\k)- 



uj 5(k) 
5*(k) u 



(5) 

where is the energy eigenvalue derived earlier. The 
density of single particle states would then be given by 



N{E) 



^Y.Jm[G^k,u)] 
4Ek^"^[G°(k,c.)] 



iiE> Ep 
iiE <Ef 



(6) 



Fig.© shows the density of states (DOS) for pure 
graphene (Fig. (jlap )and also the real and imaginary part 
of the Green's function fFig. (jlbp ). 
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(a) 



(b) X. 



FIG. 2: Feynman diagrams for the (a) Born and (b) 
T-matrix approximations for the self energy. The "x" . 
dotted line and line with arrow signify the impurity, 
impurity potential and the Green's function 
respectively. 



III. T-MATIRX ANALYSIS 




In the presence of impurities, the free particle Green's 
function gets modified. We shall, at this stage, look at 
the effect of a single impurity on the Green's function. 
We shall employ the T-matrix approach which gives ex- 
act result in the case of single impurity. As from the 
scattering theory, the modified Green's function based 
on the diagrams in Fig. ([2]) is be given by: 

G.,(u;) = Gliuj) + GUw)T{lo)G%{u) 

where the T-matrix is given as T{uj) = V/[l — VG^j{oj)], 
V is the strength of the impurity. The index /(= A,B) 
and denotes the impurity site and the lattice sites 
respectively. 

Information regarding additional poles is contained in 
the T-matrix. This implies that there are possibilities 
of existence of new states if there is an intersection of 
i?e[Gjj(w)] and 1/V. The Imaginary part is small but 
finite. If the strength of impurity is small or even moder- 
ate, there shall be no intersection with the 1/V line. It is 
only when the strength of the impurity is large, some in- 
tersection is expected and this might lead to a new state. 
Therefore, for some value of w = w*, a resonance state 
occurs at the impurity site. This new state would be lo- 
calized at the impurity site. The full Gren's function is 
given by: 



G(k,Lo) 



[G"{k,io)] 
- V[GU^) 



(7) 



Around cither K or K points (Dirac points), the G'jj 
has the following form. 



El 



""[ln(-"^'=' 



-) + in] 



(8) 



where i^k is the energy, kc is the momentum cutoff, up 
is the Fermi velocity and / is the impurity site. One 
should note that the free Green's function contributes 
equally for an impurity sitting on either sub-lattice A or 
B. Also, as, w — >■ 0; the w term in front provides a cut-off 



FIG. 3: (color online) Density of states for a graphene 
sheet with random on-site impurity potential. For, 

W=2, it is still linear. However, for W=5, it becomes 
nearly fiat. Energies are measured in the units of t. 



to the log divergence. Imaginary part gives the density 
of states as already mentioned above. Real part of the 
above expression is small. 

If we introduce a magnetic impurity at the some site 
/, there could again be appearance of new states for suf- 
ficiently large impurity strength. However, the intersec- 
tion could now be either above or below the to axis. How- 
ever, since the nature of intersection is same as for the 
single non-magnetic impurity, the physics is expected to 
remain same. 



IV. DISORDERED GRAPHENE 

We will now include diagonal (on-site) disorder and 
consider the following tight-binding Hamiltonian: 



H 



E 

<ij><j 



(9) 



on a honeycomb lattice. Here the random on-site energies 
Ci are chosen from a uniform distribution on —W/2 < 
Ei < W/2. The second term is the hopping term with 
summation over nearest-neighbour pairs of sites. From 
the matrix realization of the above Hamiltonian on a fi- 
nite lattice with periodic boundary condition in both di- 
rections, we numerically obtain the eigenfunctions and 
the eigenvalues Ei of the Hamiltonian in Eq. ^ by exact 
diagonalization of the Hamiltonian matrix. One should 
note that the spin index simply runs through the calcula- 
tion and shall be considered while making comparisions 
with the experiments. In the following, we have set the 
hopping term i = 1 as the unit of energy. 

From the exact eigenvalues we obtain the local density 
of states (LDOS) using: 



1 



E 



[E - Eu)^ + 772 



(10) 
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where N is the total number of lattice points in the sys- 
tem. 

Fig-© shows the LDOS for various disorder strengths. 
Clearly, the van Hove singularities in the pure graphene 
DOS are softened with increasing disorder strength. For 
very small disorder strength, the DOS is still linear. How- 
ever, ior W > 5, the density of states becomes indepen- 
dent of w. This could be understood on the basis of inter- 
section of Re[G^j{w)] and 1/V, as argued in the previous 
section. 



V. FINITE-FREQUENCY CONDUCTIVITY 



Wc arc interested in calculating finite-frequency con- 
ductivity for graphene. Wc use the Kubo-Greenwood 
formula4^i^ for conductivity which uses the eigenvalues 
El and the eigenfunctions (j>\ calculated earlier. The 
Kubo-Greenwood formula is given by: 



ojNuA 



where, the current matrix element AIj^^_^^{uj) is given by: 



dEj2\MLi^)\''SiE - Era)5{E + hio^ E^) 



(11) 



ieat 



lins t cosQii for n = X 



sinQij for IJ. = y 



(12) 




FIG. 4: Nearest neighbors to a Carbon atom on one of 

the sublattice on a honeycomb lattice. The dots 
represent C atoms. Qij is the angle between the X-axis 
and the other nearest neighbors. 



A. W=5 

When the density of states is independent of w (re- 
fer Fig. [21 W=5), the energy levels are nearly equally 
spaced. For a given random distribution of on-site impu- 
rity potentials, the contribution to the matrix elements 
from each energy pair of eigenvalues with fixed energy 
difference, huj = Ei — Em could simply carried out by 
keeping / — ni fixed. We average the total contribution 
by the number of such paired states considered and the 
number of random configurations, typically 5000 samples 
of the matrix elements in all. 

The expression for the finite-frequency conductivity 
then reduces to: 

cr„(w) = 8'KhNMl{uj)N {0)^/3^/31^0^ (13) 



where \l) and \m) are the eigenstates with E^ < E < 
El and O^ is the angle between the x-axis and the two 
other nearest neighbors at each site as shown in Fig. ([4]). 
Summation ,< ij > is over nearest neighbors. Here iV„ 
and A = 3V3a'^/2 are the number of unit cells and 
of each unit cell, respectively. 

The summation over Z,m in Ea. dlip is carried out by 
considering each pair of states |^) and \m) such that the 
energy difference, huj = Ei — Em, remains fixed. 

Wc shall first consider the case for = 5. 



where M^luj) is the averaged matrix element squared and 
iV(0) is the averaged density of states per site at the Dirac 
points. 

We now define the normalized conductivity 

E„(w) = 3V3hW^an{uj)/87re^ (14) 

which is convenient for plotting the data. Fig. (O is the 
data for a (20,24) lattice for various values of W. At 
high frequency the normalized conductivity is approxi- 
mately one. For higher W, the conductivity falls off at 
low frequency. This reduction of the low-frequency con- 
ductivity can be ascribed to the effects of localization. 
We show in Fig.® the normalized conductivity versus 
frequency for different system sizes at fixed W(W = 5). 
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FIG. 5: (color online) The normalized conductivity, 
E(a;) versus lo for various impurity strengths for a 

(20,24) system. In all the cases, the DOS is flat. The 
curves are drawn as a guide to the eye. Some data 
points are omitted for clarity. 
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are then renormalized to preserve the physical quantites 
as the lattice parameter is varied. We consider specif- 
ically two systems with different lattice spacing (a and 
a ) but same physical size. For these two systems to 
represent the same physical problem with different mi- 
croscopic length scales, we demand that the Hamiltonian 
parameters be so related that the physical properties are 
preserved. The appropriate physical quantites to be pre- 
served for the localization problem are the one-electron 
DOS at the Fermi energy, A^(0) and the low-frequency 
electrical conductivity^^. We define the dimensionless 
conductance which depends only on W . 

gn{uj) = fian{uj)/e^ (15) 

Then, the one-parameter RG recursion relation is: 

gn{uj,w ) = gn{uj,w) (16) 

The other MRG relation: 

N^{t ,w')^N,-,{t,W) (17) 

fixes the absolute magnitude of t to preserve the DOS. 
For weak disorder, the weak localization scaling theory 
result in 2D is^i 

91ng/ain(L) = /3(.g) = -l/27r25 (18) 

where, L is the length scale. For W = b, wc com- 
pare the (20,24) lattice with (12,14) lattice and find that 
Ag/Aln(v^) = -0.052 which matches well with the 
above value: -l/27r2 = -0.051. Thus, for W = b, the 
scaling theory of weak localization is completely obeyed. 



FIG. 6: (color online) The normalized conductivity, 
S(a;)versusaj for various system sizes. In all the cases, 
the DOS is flat (W=5). The curves are drawn as a guide 
to the eye. Some data points are omitted for clarity. 



We observe that the conductivity decreses with increas- 
ing system size which exhibits the rcnormalization group 
at work. 

We also note that, despite the fact that the lattice 
is not symmetric about x and y axes, the conductivity 
behaviour turns out to be similar in both directions (not 
shown). This is due to the fact that the spectrum of 
graphcne is conical. 

1. Macroscopic Renormalization Group 

We now turn to the macroscopic renormalization- 
group (MRG) calculation. We start with a certain set 
of Hamiltonian parameters (t and W) and calculate the 
appropriate macroscopic physical quantites as a func- 
tion of lattice parameter. The hamiltonian parameters 



B. W=2 

We now turn our attention to the case when the den- 
sity of states is not flat (e.g. for W=:2, refer Fig. [3]). The 
earlier method of obtaining the average current matrix el- 
ements squared does not work here since the LDOS is lin- 
ear. We, therefore, average the current matrix elements 
squared, |Af;mP over all I and m states such that the en- 
ergy difference fiu = Ei — Em lies in a given range (bin- 
ning). We have checked that for 14^=5 the < |M/„ip > 
obtained by this method matches with the earlier one. 

For weak disorder strength, we study the d.c. limit 
of average current matrix elements squared for different 
locations in the band. We consider a narrow band of 
states centred at different energy eigenvalues and calcu- 
late their contribution to < |Af/m(a;)p > in the limit 
a; — >■ 0. It could be directly seen from the Fig.([7{a)) that 
in the limit w — > for W=b, contribution from all the en- 
ergy states to conductivity is nearly equal. This implies 
that all the states are of the similar nature. However, 
in the same limit for W=2, the states close to the Fermi 
energy have a larger contribution and the contribution 
decreases as we move away from the band centre as seen 
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FIG. 7: (color online) Energy resolved current matrix 
element squared for(a) W=5 and (b) W=2 for a (16,20) 
graphene system subject to disorder. 



in Fig. (HKb)). This implies that, close to band centre, 
the states are more extended (compared to the states 
away from the band centre). We shall return to the issue 
of localization in the next section. 

We may now ask : what is dependence of | Mim{u! — >■ 
0) p for W=2 on energy eigenvalues as we move along 
the band? We polynomial fit the data and find (Fig.®) 
that the behaviour is inversly proportional to the square 
of the energy eigenvalue. The dependence is of the form 
a/Ef + b. For W =2 and system size (16,20), we find 
that a = 0.0010, and 6 = 0.0021. 

We study the scaling behaviour of the constants aN 
and bN with increasing system size. (Fig. fTU]) . We have 
plotted the behaviour of aN and bN versus 1/N. It could 
be seen that there is not a significant change in the value 
of aN and bN as we increase the system size and the 
dependence of these quantities on 1/N is almost linear. 
As TV 0, the values of aN and bN are 1.06157 and 
3.8294 respectively. 

The first term (~ ^/Ef) is what is expected from the 
Boltzmann theory. The other term, independent of Ei 
is the non-Boltzmann contribution. A particle with the 
Hamiltonian described by Eq.([T]), which supports both 
positive and negative energy states (positive and negative 
bands), could be in a linear combination of the positive 
and negative energy eigenstates. The Boltzmann term 
arises out of transition between term in the same band. 
Off-diagonal terms in the velocity operator correspond 
to the transition between the same k but belonging to 
different bands and is responsible for the non-Boltzmann 
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FIG. 8: (color online) Energy resolved current matrix 
element squared for W=2 for (a) (16,20), (b) (20,24) 
and (c) (24,30) graphene systems subject to disorder. 



contributions to the conductivity. The origin of this term 
could also be understood from the fact that the veloc- 
ity operator is non-diagonal in the helicity basis. An 
analytical calculation for Boltzmann conductivity with 
screened Coulomb scatterers and incorporating the off- 
diagonal terms in the velocity operator was carried out 
by Trushin et al^ 



1. Conductivity and Mobility 

In order to study d.c. conductivity, we go back to 
Ea. (fTT|) whereby d.c. conductivity could be written as 
the product of the current matrix elements squared in the 
limit w — )- and DOS squared. In order to do this, we 
polynomial fit the DOS squared and the average current 
matrix elements squared in the d.c. limit. 

For W =2, we have seen that DOS ~' lo. Carrier con- 
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FIG. 9: (color online) The matrix element squared, 
\M{uj — )■ 0)1 v/s El for a (16,20) graphene lattice and 
W=2. The curve is a a/Ef + 6 fit with a= 0.0010 and b 
= 0.002 
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FIG. 10: (color onhne) Scaling of aN and bN with 
respect to system size. 



centration, n is defined as n = 3.8 x 10^** Ep{cm~'^). 
Putting in the values of the parameters, t — 2.8eV and a 
= 1.42 A, together with the the bahavior of | A/;m(aj — > 
0) 1^ versus Ei (obtained in the previous subsection) and 
n versus Ei , we obtain conductivity as a function of the 
carrier concentration. Here the factor of two for spin mul- 
tiplicity has already been taken into account. For large 
carrier concentration, wc notice that as Ei increases, the 
constant term ( " 6" ) in the fit dominates the current ma- 
trix element contribution. This together with the DOS 
shows a linear dependence on the carrier concentration. 

We also note that for n a; zero, the conductivity does 
not go to zero. It turns out that for small W, both the 
|M;mp and the DOS together conspire to give a non- 
zero d.c. conductivity (cTmin), independent of Ep, close 



FIG. 11: (color online) Behavior of conductivity with 
carrier concentration for W = 2. The minimum value of 
conductivity at n = 0, (Tmin is about 1.3e^//i. 



to the Dirac point. Comparision of different system size 
suggests that as the system size is increased the lesser 
number of cigcnsatatcs contribute to (Tmin for very small 
doping concentration. However, the ~ ^/Ef dependence 
is still retained. This translates to the fact that as we in- 
crease the system size, the 1 / Ef shifts towards the Fermi 
energy as seen in Fig.®. 

The value of this minimal conductivity, (Jmin ~ 4e^/7r/i 
in confirmity with the previous results obtained with 
SCBAii^i^. 

We also compare this behaviour with respect to the 
system size and find that for larger system size, the num- 
ber of such extended states decreases as shown in Fig. ([8]) 
but still retaining its ^ ^/Ef dependence. 

Thus, away from the zero doping concentration, one 
obtains a linear dependence of conductivity on the doping 
concentration in confirmity with the experimental find- 
ings. However, there is a quantitative difference between 
the two. The slope of the a v/s n (Fig.(|ll|)) is very small 
compared to the experimental results. 

The dependence if a.c. conductivity on w can be ob- 
tained from the energy resolved matrix elements squared. 
For a fixed energy, the dependence of |M;m(a;)p ^ l/w^. 
This together with the 

Mobility is defined as: fi = a/ne. Since, a goes linearly 
as n, n should nearly be a constant for energies away from 
the Dirac points. Since n oc Ep, from Ea. (fT3)) . we find 
that fi cx |M;m(a;)p. 



3V3/1 X 3.8 X 10i4^a;2 



b) 



(19) 



In Fig. (fT2|) . we plot versus gate voltage, Vg where n 



7.2 X 10 V^. We find that the mobility rises sharply close 
to the Dirac point. This sharp rise around the Dirac 
points is in a very good qualitative agreement with the 
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FIG. 12: (color online) Plot of mobility v/s gate voltage 
(carrier concentration) for a (20,24) graphene lattice for 
different disorder strengths W. The plot matches 
qualitatively with the experimental result. 
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FIG. 13: (color online) Average current matrix elements 
square v/s W for different system sizes. The y-axis 
corressponds to the average dc value of for a band of 
ten states around the fixed energy 0.2t. 



is not sensitive enough to capture signatures of localiza- 
tion for small W. Therefore, we, again, turn towards the 
energy resolved studies of averaged current matrix ele- 
ments squared, (Fig. [7|). As mentioned earlier, for a few 
states close to the band centre have a very large averaged 
current matrix clement squared value. This implies that 
these state are extended. We notice that this is the case 
independent of the system size. Also, the experimental 
finding of a universal zero-bias minimum conductivity, 
fmin suggests that the states close to the band centre 
should be extended, as we find which is in conformity 
with the findings of Amini et al^. This, however, is in 
contradiction with the findings of Xiong et al^. 

One might, also, question the | Mim{oj) ^ dependence 
on the energy eigenvalues, Ei since it could very well be 
a finite-size effect. For low disorder strength, the local- 
ization length, ^ varies exponentially with the inverse of 
disorder strength i.e. ^ exp{l/W) and in order to see 
the effects of localization, one needs to change the system 
size exponentially. We note that the ^ ^/Ef behaviour is 
present at all system sizes starting from system contain- 
ing nearly 300 lattice points((8,10)) to nearly 3000 lattice 
points((24,30)). 

Also, that for W=3, the |Mim(a;)p is independent of Ei 
(not shown) suggesting a disorder induced crossover from 
localized to delocalized phase for 2 < W < 3. We check 
this by studying the averaged current matrix element 
squared over a range of W for different system sizes, at a 
fixed energy. We have fixed the energy at 0.2t and have 
considered ten states around this energy. We study the 
contribution of these ten states towards the \Mim{i^)\'^ ■ 
The results are shown in Fig. lfTS]) . The inset of Fig. (fT3|) 
shows the plot between averaged current matrix element 
squared over a range of W for three different system sizes 
(16,20), (20,24) and (24,30). The cross-over is highlighted 
in the main panel which shows the logarithm of averaged 
current matrix elements squared versus W. We indeed 
find a crossover around W^=3. However, this crossover is 
not between localized and delocalized phases but from 
localized to a very weakly localized behaviour. This is in 
agreement with the scaling theory of localization which 
predicts that all states in 2D are localized. 



experiments. However, the values of mobility is around 
1000 times lower than the experimental results. 

We would like to add that, in the adopted formalism, 
it is impossible to distinguish the effect of inter-valley 
scattering or to determine it's contribution towards con- 
ductivity. However, since we are dealing with very short 
range scatterers, it would be impossible to avoid large 
momentum interactions and hence intcr-valley scatter- 
ing. 

VI. LOCALIZATION 

In order to study localization, we looked at IPR ver- 
sus energy eigenvaluei(not shown) and found that IPR 



VII. COMPARISION WITH SQUARE LATTICE 

We, now, compare the normalized conductivity, S(w) 
between square lattice and graphene ioT W — 5. This 
comparision has been shown in Fig. lfn)) . We observe that 
the conductivity matches well at large frequencies for 
both the lattices. However, at low frequencies, the gra- 
hene has a lower conductivity suggesting that the states 
are more localized. This could mean that the states in 
graphene are more susceptible to localized when subject 
to disorder, as also suggested by Xiong et al^ 

We also compare the intensity plots for graphene and 
square lattice for different values of disorder strengths 
and fixed system size. The results are shown in FIG (|15|) . 
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FIG. 14: (color online) Comparision of |Af p values 
between a square lattice of size (36,36) (=1296 lattice 
points) and graphene of size (16,20) (=1280 lattice 
points) for W=5. Averaging has been done over 50 
configurations of disorder. 



We find that there is no signature fo puddle formation 
for square lattice. The reason for the formation of puddle 
for low carrier concentration in graphene is yet to be 
understood. 



VIII. CONCLUSION 

In conclusion, we have studied the finite-frequency 
electrical conductivity in graphene under diagonal (on- 
site) disorder of various strength and in the presence of 
inter-valley scattering using Kubo-greenwood foumula. 
For moderate disorder strength, we find that for differnt 
system sizes, fixed W, scaling theory is at work. We 
made contact with the weak localization result of Abra- 
hams et alM- and found a very good agreement which 
means that for W = 5, logarithmic scaling is obeyed. 
We compare normalized conductivity of graphene with 
that of a square lattice and find that graphene is more 
susceptible to localization when subject to disorder. We 
have established that, for low disorder strength, W = 2, 
the states away from the band centre are more localized 
comared to the ones close to the band centre. For low 
disorder strength, we have calculated the conductivity for 
low-disordered graphene and have found the results are in 
disagreement with that of Boltzmann conductivity. Also, 



the conductivity and mobility are in qualitative agree- 
ment with experiments. Also, for weak disorder strength, 
it is the competition between |M(a;)p and DOS which 
gives rise to a universal conductivity minimum. The lin- 
ear rise of conductivity with carrier concentration is due 
to a term unaccounted for in the Boltzmann expression 
for conductivity. We have also established that, for low 
disorder strength, W = 2, the states close to the band 
centre are extended and that there exists a crossover at 
W ^ 3. Comparative study of intensity plots for states 
close to the band centre for graphene with square lattice 




FIG. 15: (color online) Intensity plot of energy 
eigenstate next to the band centre for pure graphene 
sheet ((a),(c) and (e)) and square lattice ((b), (d) and 
(f))for different values of disorder strength viz. 
disorder- free (pure) ((a) and (b)), W=2 ((c) and 
(d))and W=5 ((e) and (f)). The system sizes considered 
are comparable for both the systems, consisting of 8000 
((40,50)) and 8100 ((90,90)) lattice points. 



shows clear signatures of puddle formation in graphene. 
Although we have studied a simple disorder model in 
graphene, some of the feature studied could be generic. 
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